Mott insulators in an optical lattice with high filling factors. 



(N 
O 

o 

(N 

» 

o 
O 
in 

(N 



CO 



I 

C 

O 

o 



> 
O 

in 
o 

<N 
O 

c3 



-a 
c 

o 
o 



X 
S3 



D. van Oosten, 1 - 2 P. van der Straten, 2 and H. T. C. Stoof 1 

1 Institute for Theoretical Physics, Utrecht University, 
Leuvenlaan 4, 3584 CE Utrecht, The Netherlands 
2 Debye Institute, Utrecht University, Princetonplein 5, 3584 CC Utrecht, The Netherlands 

We discuss the superfluid to Mott insulator transition of an atomic Bose gas in an optical lat- 
tice with high filling factors. We show that also in this multi-band situation, the long-wavelength 
physics is described by a single-band Bose-Hubbard model. We determine the many-body renor- 
malization of the tunneling and interaction parameters in the effective Bose-Hubbard Hamiltonian, 
and consider the resulting model at nonzero temperatures. We show that in particular for a one or 
two-dimensional optical lattice, the Mott insulator phase is more difficult to realize than anticipated 
previously. 
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Introduction. — The behaviour of trapped Bose- 
Einstein condensates offers a large amount of interesting 
features. Specifically, the phase coherence of a conden- 
sate creates the prospect of various interference experi- 
ments, as shown for the first time in an experiment per- 
formed by Andrews et al. fjj]., In another experiment a 
large number of condensates trapped in a periodic lattice 
potential have been made to interfere 0. Having shown 
experimentally that in general condensates are phase co- 
herent and thus show off-diagonal long-range order, the 
question arises if this long-range order can also be de- 
stroyed in a controllable way. This was only very recently 
achieved in a beautiful experiment by Greiner et al. || . 

In this last experiment, a trapped Bose-Einstein con- 
densate is put into a three-dimensional optical lattice. 
The number of atoms in the condensate is sufficiently 
large to obtain a filling factor of almost two atoms per 
site. By increasing the intensity of the lattice light, a 
quantum phase transition from a superfluid state to a 
Mott insulating state is achieved. In the insulating phase 
all phase coherence is lost due to quantum fluctuations. 
The transition was predicted to occur in this system by 
Jaksch et al. (J] , and the observed critical conditions for 
the transition are in good agreement with the results of 
an on-site mean-field theory H. This indicates that the 
single-band Bose-Hubbard model used in Refs. |j| can 
accurately describe a gas of ultracold bosonic atoms in 
an optical lattice when the filling factor of the lattice is 
of the order of one. 

However, this model is no longer valid in the case of 
higher filling factors such as described in the experiments 
°f H 01- The theories mentioned above use single-particle 
wave functions corresponding to the lowest band of the 
lattice to calculate the microscopic parameters of the 
single-band Bose-Hubbard model as a function of the lat- 
tice parameters. In the case of high filling factors more 
than one band is generally populated, leading to a multi- 
band Bose-Hubbard model. The interaction effects that 
occur under these circumstances have not been consid- 
ered previously. Furthermore, the effects of thermal fluc- 



tuations are also not understood, even in the single-band 
Bose-Hubbard model. There are studies that describe 
number squeezing in an optical lattice at nonzero temper- 
ature J|, ||] , but they are not in the strongly-interacting 
limit that is relevant for the Mott insulator. Our main 
objective here is to develop an effective theory that can 
deal with these issues. 

Solving the multi-band Bose-Hubbard model. - To 
solve the problems associated with high filling factors, 
we have to deal with the many-body physics at every 
site. Since the high filling factors of interest are experi- 
mentally most relevant in low-dimensional lattices ^ Q , 
we discuss the energy scales involved in those systems. In 
a low-dimensional lattice, we can approximate the on-site 
trapping potential by an anisotropic harmonic potential. 
The oscillator frequencies uu and ui±_ correspond to the 
trapping frequencies in the directions parallel and per- 
pendicular to the periodicity of the lattice, respectively. 
Because the typical size of a lattice well in the paral- 
lel direction(s) is much smaller than in the perpendicu- 
lar direction(s), we immediately have that hu>± <C Hlou. 
Furthermore, for the experimental conditions of inter- 
est ||, the temperature is in between the two trap- 
ping frequencies, i.e., huj± <C ksT <C tiu)\\. This implies 
that in every site the gas is in the parallel direction(s) 
in the ground state of the potential, but that it occu- 
pies many states in the perpendicular direction(s). In 
particular, this holds for the thermal cloud of the gas. 
As a consequence, the effective dimensionality of the gas 
at every site is reduced and the thermal excitations are 
only present in the perpendicular direction(s). It is un- 
der these conditions that we are able to solve the relevant 
multi-band Bose-Hubbard model by using the following 
two-step procedure. 

First, we solve the many-body physics at every site. 
Due to the famous infrared problems of a one or two- 
dimensional Bose gas this is not an easy task, and an 
accurate equation of state for these gases was developed 
only very recently in the weakly-interacting limit [ pd| . 
This equation of state is found by treating phase flue- 
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tuations in the (quasi)condensate exactly and we can in 
particular use it to determine at every temperature the 
number of atoms in the (quasi)condensate N (T). Fur- 
thermore, it is shown in Ref. JTT| ] that even in the pres- 
ence of phase fluctuations, the Gross-Pitaevskii equation 
can still be used to calculate the density profile of the 
(quasi)condensate. 

Secondly, we consider the coupling between the sites. 
Since we have a (quasi)condensate at every site, the 
coupling between sites will be dominated by tunneling 
from (quasi)condensate to (quasi)condensate as opposed 
to (quasi)condensate to thermal cloud. This means that 
we can describe the coupling between sites by a single- 
band Bose-Hubbard model. The important parameters 
in the Bose-Hubbard model are the on-site interaction 
energy U and the energy t associated with the tun- 
neling of atoms between nearest-neighbor sites. Both 
energies can be calculated from the knowledge of the 
(quasi)condensate wave function Vo( x ) = ^/ n o( x )e 4,? , 
where no(x) is the density profile and d is the global 
phase of the (quasi)condensate. The interaction energy 
U is proportional to J dx|-0o(x)| 4 , whereas the tunneling 
energy t requires the evaluation of an overlap integral be- 
tween the (quasi)condensate wave functions of two neigh- 
boring sites in the parallel direction(s). 

To describe the effect of the interatomic interaction, 
we thus need to determine how the (quasi)condensate 
wave function changes as a result of the on-site inter- 
actions. Since the mean-field interaction obeys Tiuj±_ <C 
NqU -C TktMi under the experimental conditions of in- 
terest, we can write the three-dimensional wave func- 
tion of the condensate as a product of a single-particle 
ground-state wave function in the parallel direction(s) 
and the (quasi)condensate wave function in the perpen- 
dicular direction(s). If we substitute this product wave 
function into the Gross-Pitaevskii equation and integrate 
out the parallel direction(s), we arrive at an effective 
equation for the (quasi)condensate wave function. Be- 
cause of the above mentioned inequality, we can subse- 
quently solve this equation using the Thomas-Fermi or 
local-density approximation jl2) . To quantify the differ- 
ences between the (quasi)condensate wave function and 
the single-particle ground-state wave function, we define 
a dimensionless coupling constant g both in the noninter- 
acting and in the interacting case. The first parameter 
we call the bare coupling constant <?b = E/bAb and it 
is calculated with the single-particle ground-state wave 
function in every site. The second parameter we call the 
renormalized coupling constant <?r — £/ R /i R and it is 
calculated using the single-particle ground state in the 
parallel direction(s) and a Thomas-Fermi density pro- 
file in the perpendicular direction(s). Because we have 
already included the on-site interaction effects in this 
coupling constant, we can now write down a renormal- 
ized single-band Bose-Hubbard model for the total opti- 
cal lattice, where the creation and annihilation operators 



a], and hi, respectively, and the number operator hi are 
not associated with the Wannier states of atoms in the 
lattice, but with the macroscopic wave function of the 
(quasi)condensate in each site. In particular, we have 

H = -t R af Qj + hj (hj - 1) - /iR n^l) 

(id) 

where /xr is the effective chemical potential. The inter- 
action parameter is given by £/r = d 2 F os /dN 2 \ N _ N = 
9fi os /dN\ N _ N , where F os is the on-site free energy and 
/i os is the on-site chemical potential. Formally, the effec- 
tive chemical potential is given by /xr = /i — /i os — E/r/2, 
where the last term is substracted from the interaction 
energy to make the analogy to the single-band Bose- 
Hubbard model of Refs. |5) complete. 

It is important to understand that the hopping term 
only describes hopping between the (quasi)condensates 
in neighbouring sites. While it is clear that this is a very 
good approximation in the case of neighbouring conden- 
sates, it may not be immediately obvious in the case of 
neighbouring quasicondensates. However, it should be 
noted that the effect of the hopping is only large when 
the system is in the superfluid phase, in which case all 
the sites couple to form a true three-dimensional con- 
densate. The tunneling strength can be calculated in the 
tight-binding limit and depends only on the overlap in 
the parallel direction(s). As a result the bare and renor- 
malized values of t are equal. However, the interaction 
energy is strongly reduced due to the repulsive on-site in- 
teractions which spread-out the (quasi)condensate wave 
function considerably. We find for a two-dimensional gas 
that 




and for a one-dimensional gas that 




Here a is the positive s-wave scattering length of the 
atoms, In = y/h/moJu and £± = y^h/muj± are the har- 
monic oscillator lengths in the parallel and perpendicu- 
lar directions, respectively, and Rtf is the Thomas- Fermi 
radius of the (quasi)condensate. The physical interpreta- 
tion of Eqs. (||) and (||) is that as a result of the repulsive 
interatomic interactions, the (quasi)condensate reduces 
its total energy by increasing its size in the perpendicular 
direction(s). This can be seen from the fact that the de- 
crease of the coupling constant is inversely proportional 
to the increase in the surface or length of the two or one- 
dimensional gas, respectively. Note that his reduction is 
particularly important for the critical conditions, which 
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can be written as i?r > 4zN for large No || . In order to 
verify the consistency of our two-step approach, we ex- 
plicitly check the relevant energy scales using our results 
of Eqs.(^) and (Q). First, we assumed that the mean- 
field interaction energy NqUh is much smaller than the 
trapping frequency in the parallel direction(s) 7ko\\. This 
requires that for a one-dimensional lattice 
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(4) 



Second, we also assumed that the crossover temperature 
for the formation of a (quasi)condensate in two dimen- 
sions is much lower than tkuu. This results in 



(5) 



where N is the total number of atoms at every site. 
For typical numbers used in the experiments by Orzel et 
al. ||, we find the condition Nq < N -C 10 5 , which means 
that our assumptions are valid for even the largest fill- 
ing factor reported. For the case of the two-dimensional 
lattice of Greiner et al. 0, a similar inequality can be 
derived. It is found that this experiment is also in the 
regime where our assumptions are valid. Note that the 
use of the Thomas-Fermi approximation also imposes a 
lower limit on the filling factor, namely, N a 3> l±. How- 
ever, when the filling factor is below this limit, we are in 
the regime where we can safely use the bare coupling 
constant. We thus conclude that depending on the fill- 
ing factor, either our renormalized or the bare theory is 
applicable to these experiments. 

Thermal effects. - Besides the effect of thermal 
fluctuations on the number of (quasi)condensate atoms 
per site N (T), which is accounted for by the equa- 
tion of state of the low-dimensional Bose gas [jlT) , there 
is also the effect of thermal fluctuations on the renor- 
malized single-band Bose-Hubbard model itself. These 
thermal fluctuations are also present in a lattice with 
low filling factor. To study these thermal fluctuations, 
we use an on-site mean-field Hamiltonian that we can 
derive from Eq. ([!]), using the approach presented in 
Ref. ||. We decouple the tunneling term, by intro- 
ducing a complex mean field parameter ip as follows 
a\&j — ipa,j + a\ip* — \ip\ 2 ■ Physically, ip is the super- 
fluid order parameter which we choose to be real in the 
following. Performing the above substitution, wc find 



H{ip) — ~zti)j{a) + a) H — — n(n — I) — /iR,n + ztip 2 , (6) 

with z the coordination number. Since this is an 
on-site Hamiltonian, we have dropped the site indices 
for simplicity. Moreover, we assumed that the chemi- 
cal potential is chosen such that the expectation value 
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FIG. 1: Phase diagram of the Bose Hubbard model in terms 
of the chemical potential nn/zt and the coupling constant 
Un/zt. The solid and dotted lines correspond to T = and 
T = O.lUc (where U c is the critical Ur for the iV = 100 lobe), 
respectively. The inset shows a qualitative phase diagram in 
terms of the temperature T and the coupling constant g. N, 
SF and MI indicate the normal gas phase, the superfluid and 
the Mott insulating phase, respectively. 



of the number operator n is equal to the number of 
quasi(condensate) particles N (T) in every site. 

The zero-temperature phase diagram of this mean-field 
theory can be solved exactly |], [lOj and is shown in Fig. [jj, 
where the Mott insulator phases correspond to the var- 
ious lobes. For nonzero temperatures the model can no 
longer be solved analytically and we have to resort to nu- 
merical methods. If we put ■0=0, we find that the eigen- 
states of Eq. (|]) are given by number states. Using a basis 
consisting of these number states, running from a certain 
minimum filling factor iV m i n to a certain maximum N max , 
we can calculate the grand-canonical partition function 
Z(tp) = Tr[e- H W/ fcflT ] by diagonalizing the mean-field 
Hamiltonian given in Eq. (^). Next we determine the 
thermodynamic potential £l(ip) = — k^T hi Z(ip) as a 
function of the order parameter if. For zero tempera- 
ture the calculation converges when iV max ~ Nq + 4 and 
iVjnin ~ iVo — 4, where N$ is the filling factor of the 
relevant Mott-insulator lobe. For nonzero temperatures, 
more states must be included. 

To obtain the relevant thermodynamic quantities, we 
minimize the grand potential fi(V0 and the value of tp 
at the minimum of corresponds physically to the 

square root of the number of atoms that is superfluid 
in the direction(s) parallel to the periodicity of the lat- 
tice. In the Mott insulator the gas is only insulating in 
the direction(s) parallel to the periodicity of the optical 
lattice, whereas it is always a superfluid in the perpen- 
dicular direction(s) . The other quantity of interest is the 
value of the number fluctuations. This number is im- 
portant because in the Mott-insulator phase the number 
fluctuations are exactly zero. 

The final results of the calculations are shown in Fig. |^. 
In these figures, the longitudional superfluid fraction and 
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FIG. 2: Condensate fraction (a) and particle number fluctu- 
ations (b) as a function of the coupling constant, for fcsT = 
0,0.08?/ and 0.145C7 (solid, dashed and dotted lines respec- 
tively). The coupling constant Un/zt = 440 (the dashed line 
in Fig 0). 



the number fluctuations are plotted along the dashed line 
in Fig. [I] for different temperatures. It can clearly be 
seen from Fig. |2|(a) that the superfluid part of the phase 
diagram decreases with increasing temperature. In addi- 
tion Fig. ||(b) shows that at zero temperature the den- 
sity fluctuations drop exactly to zero in the Mott insu- 
lating regions, but that this does not happen at nonzero 
temperature. This is a result of the fact that the su- 
perfluid to Mott insulator transition is a quantum phase 
transition. The reason that there is still a reduction in 
particle-number fluctuations at nonzero temperature is 
that the excitation spectrum of a fluctuation is gapped 
in this region ||, which means that the fluctuations are 
exponentially suppressed. Due to this strong suppression 
of the number fluctuations, one will be able to observe 
a phase which is formally not a Mott insulator, but ex- 
perimentally has very similar features. Another feature 
we can clearly see in Fig. ||(b) is that part of the phase 
diagram where the number fluctuations are suppressed 
also decreases with increasing temperature, and shrinks 
in the opposite direction of that of the superfluid part. 

On the basis of the above calculations, we can draw 
the nonzero temperature phase diagram shown in Fig. [I]. 
In this figure, the solid lines indicated the superfluid to 



Mott insulator transition at zero temperature and the 
dotted lines indicate the superfluid to normal transition 
at nonzero temperature. The inset shows the phase dia- 
gram in terms of the temperature and the coupling con- 
stant. This diagram agrees very well with the general 
description given by Sachdev (h}. 

Conclusion - We have shown that for low-dimensional 
lattices, which generally have a filling factor much larger 
than one, we should in principle solve a many-band Bose- 
Hubbard model. This can be achieved by first solving the 
on-site many-body problem, and then deriving an effec- 
tive theory that describes the coupling between the sites 
in the optical lattice in terms of a renormalized single- 
band Hubbard model. We have calculated the effects of 
thermal excitations in this renormalized model and we 
have shown that the number fluctuations in the above 
model can only drop to zero in the absence of ther- 
mal fluctuations. However, if the temperature is suffi- 
ciently low, the number fluctuations are exponentially 
suppressed. This means that at a certain nonzero tem- 
perature, the crossover to the Mott insulator phase can 
still be observed if the coupling constant is increased to a 
value larger than the zero-temperature critical value (cf. 
Fig. [I]) . It is important to realize that to experimentally 
obtain the Mott insulator with a large filling factor No, 
the coupling constant or = U^it must be larger than 
AzNo. However, Eqs. (g) and (Q) show that the renor- 
malized coupling constant is much smaller than the bare 
coupling constant for a low-dimensional optical lattice. 
We therefore conclude that in terms of the bare coupling 
constant, which is the experimentally relevant control pa- 
rameter, the Mott insulator phase is much more difficult 
to obtain than is naively anticipated on the basis of a 
purely single-band Bose-Hubbard model. 
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